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Q\ , Hubbard model for the charge carriers in Cu02 plane. For this purpose the 



The t-J model effective hopping integral is determined from the three-band 



values of the superexchange constant J and the charge-transfer gap E gap are 

calculated in the framework of the three-band model. Fitting values of J and 

Egap to the experimental data allows to narrow the uncertainty region of the 
O 

three-band model parameters. As a result, the t/J ratio of the t-J model is 

fixed in the range 2.4 -j- 2.7 for holes and 2.5 -r 3.0 for electrons. Formation 

of the Frenkel exciton is justified and the main features of the charge-transfer 

spectrum are correctly described in the framework of this approach. 
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I. INTRODUCTION 

A large amount of works dedicated to high-T c superconductors agree that an appropriate 
electronic model which contains all essential orbitals is the three-band Hubbard model |l]-[3| . 
Some other works developing a multiband approach j|-|7J also support this belief. 

It is also widely accepted that the low-energy physics of insulating compounds can be 
described by the Heisenberg Hamiltonian. According to the earlier work by Zhang and 
Rice H and latter studies 0- |HJ lightly doped systems are described by the simple t-J 
model. Different techniques provide different exactness of the three-band model to t-J model 



mapping and generate a wide type of generalizations jTl| , [14|,[15 



In our previous works [|l2|,|r|,|T5|] a consistent low-energy reduction of the three-band 



model to the generalized t-J model in the realistic range of parameters has been performed. 



It has been shown |L5[ that the second-order corrections to the local energy of the carrier 
and its hopping integral are small (< 5%). The role of the next-nearest-neighbor terms has 
also been discussed. 

It is commonly believed that the 90% or even higher accuracy of the t-J model as the 
low-energy electronic model for high-T c superconductors justifies its wider study [|T6|-|i"8|| and 
remains this model the main pretender in describing of the superconductivity in cuprates 
P"5| [?T| . Recent angle-resolved photoemission experiments [F2^] can be interpreted as a direct 
support of some t-J model properties |23| . Some other anomalous behaviors of systems in 



normal state as well as superconductivity itself seem to be described in the framework of 
the t-J model |JH1I, @. 



The reasonable question from this point of view is: what is the role of either the three- 
band or more complex first-principles models? There are several answers: (i) calculating 
parameters for the t-J model for real systems; (ii) giving insight into the experiments in- 
cluding not only the simple t-J model degrees of freedom. 

In this paper we mainly address the question of the t-J model parameters. Superexchange 



constant J for the t-J model is directly measured p3j25|. Hence, the parameter to be 



determined is the effective nearest-neighbor hopping integral. The problem of its calculation 
is not connected with the accuracy of low-energy mapping from the three-band model, which 
is always very high, but results from uncertainty in the three-band model parameters. The 
three-band model in conventional formulation contains as inner parameters two on-site and 
one inter-site Coulomb repulsions, two hopping integrals, and split of the levels JIJ which 
are not directly measurable. Some of them are fairly bad determined. This makes the 
calculation of the hopping integral for real systems questionable and even controversial. 

We develop an obvious idea of fixing the three-band model parameters by using exper- 
imental data. This idea has already been exploited in the cluster calculations for spectro- 



scopic data and for the superexchange J in our previous work ||15| . Now, on the basis of 
better understanding of the low-energy model of the electronic system and magnetic polaron 
nature of the t-J model carriers |i8|j26| , |27[| , we calculate quite accurately the charge-transfer 
gap. Selfconsistent calculation provides a narrow range of possible values of t/J. Also, 
an excitonic feature of the charge-transfer spectrum is obtained in agreement with resent 
experiments |]2S|j . 

Paper organized as follows. In Sec. II we discuss the low-energy limit of the three- 
band model and experimentally observable quantities. In Sec. Ill the calculation of the 
superexchange constant J and the charge-transfer gap E gap is produced. In Sec. IV we 
discuss the dependence J and E gap on the parameters of the three-band model and determine 
the range of the hopping parameter t for electrons and holes. In Sec. V the properties of 
the excitonic state is considered. Sec. VI presents our results and discussions. The technical 
aspects of work are given in Appendix. 

II. THE LOW-ENERGY LIMIT OF THE THREE-BAND MODEL AND 

OBSERVABLE QUANTITIES 

Previously, it was suggested that the three-band Hubbard model is an appropriate start- 
ing point for describing the electronic structure of Cu0 2 planes 0,@]. The Cu d x 2_ y 2 orbital 



and pa(x,y) orbitals are strongly hybridizated. These orbitals are explicitly treated in the 
three-band model with the justifiable assumption that other orbitals does not directly par- 
ticipate in the low-energy dynamics. The full Hamiltonian of the model is defined by [|TJ 

H = H + H t + AH 

H = e d J2 < + e P E nP ma + u dYl n ?Xl > 

l,a rn,a I 

H t = t pd E ( d taPma + H.C.) 
<lm>,a 

~ Up E (PmaPm'a + H.C.) , (1) 

<mm' > ,a 

AH = U P Y,n p m1 n p mi + V pd £ n d la n%, (2) 

m <lm>,a/3 

in standard notation of holes at 0(p) and Cu(rf) sites. The sign convention for oxygen 
orbitals in H t |15| , p9[ is accepted. Our approach to description of low-energy properties of 



the above model []T3J is based on taking into account the main Coulomb (JJd) interaction 
exactly and the others as perturbations. 

In order to justificate this method we briefly describe here the magnitudes of the three- 
band model's parameters. Different experimental [|BH7|.|30| , atomic [31|, and band calculations 



J33| show that: U d = 5^7 eV [0, 7-5-11 eV |,|], U p = 3^8 eV §, V pd = 0-5-1.7 eV. U p 
is always less than U d . There is a general agreement for the Cu-0 system that A = e p — e d 



is always > and < U d [ 30| . It reflects the facts that the first hole in the unit cell is 
predominantly at the Cu site and the added hole has an oxygen's character, tpd — 1 -5- 1.6 
eV (and it is unlikely that it is less than 1 eV) , t pp = 0.5 -5- 0.7 eV fl, |32| , |33[| . This set of 
magnitudes will be called hereafter the realistic region of parameters. 

The consistent low-energy reduction of the three-band model to the generalized t- J model 



has been performed in previous works p~2|JT~3| , p~5|1 . Our method of the low-energy reduction 



has been based on construction of a set of local states with different number of holes over 
the filled atomic orbitals. The most essential states are: 

1. The vacuum state or the vacancy which is simply 

|u>=|0>, (3) 



2. The one-hole states which represent the ground state of Cu02 plane 

|/a >= |a >= U\da > -V\pa >, (4) 

where [dot > and |ga > are the copper and symmetrical oxygen hole states with spin 
projection a, respectively 

3. The two- hole states which are the Zhang- Rice singlets: 

\ c > = c/i|rf T d l> +^i|gt q i> 

+ W 1 (\d]qi>-\diq]>)/V2. (5) 



The coefficients U, V, U\, V\Wi are functions of the parameters of the three-band model |15 
At half filling Hamiltonian (|l]||) is reduced to the Heisenberg Hamiltonian with spin-1/2 
which are antiferromagnetically ordered due to the second-order virtual transitions through 
the set of two-hole states. Note, that the above named spins-1/2 are exactly states |/a > 
(Eq. 



It has been shown [I3,I5| that for the case near to half filling the Hamiltonian of the 
three-band model is reduced to the Hamiltonian of the t-J like model of singlets, vacancies 
and spins: 

H+.j = (E v - /i) £ XT + (E c + /x) E XT 

I I 

+ t e J2 xpx? v + th E X T*? C 

<W>,a <W>,a 

+ JE S ^'- ( 6 ) 

<w> 

Where Xf 6 = \al >< lb\ are the Hubbard operators at the site /, S; = <j a pX\ a /2. The 
constants E v and E c are the local energies of the vacancy and singlet; \x is the chemical 
potential; t e and th are the hopping integrals for the vacancy and singlet (electron and 
hole), respectively; J is the exchange constant. All five parameters E v ,E c ,t e ,th and J are 
functions of the three-band model parameters. It has been shown that relative magnitudes 



of the omitted terms in the Hamiltonian H t -,j (Eq. (||)) are of the order of ten percent 



We believe that the Hamiltonian (|6|) describes many important properties of the cuprates. 
Therefore, the real values of its parameters are of great interest. As was noted above, the 
parameters of the primary model (|]],@) are known with low precision. In this situation 
calculation of the observable quantities is an urgent issue since it provides a way to fix 
parameters of Hamiltonians (0,0) and (P). The best-defined experimental values which we 
can describe quite accurately are the superexchange constant J and the charge-transfer gap 
Eg a p. The experimental values of J are 0.14 eV and 0.17 eV for the lantan and yttrium 
systems, respectively p4|j25| . These values of J follow from measurement of the velocity of 
sound for magnons. The value of the charge-transfer gap is known from a variety of optical 



measurements p8| , |30|j33| , |34]| and is near to 2.0 eV. Observation of photoconductivity at the 
same energies shows that the excitations result in separated electrons and holes [[2l],[34|]. We 
have taken most of clear experimental features of the charge-transfer spectrum from Ref. 



2"5] where photoconductivity as well as reflectivity data for La2Cu04 are presented. 



III. CALCULATION OF THE OBSERVABLE QUANTITIES 

The expression for the AF coupling constant J in the framework of our approach is 

j = -2h 1 v% + Y,xn%r- (7) 

The first term in Eq. (^) represents the exchange energy between two holes (spins) due to 
the repulsion at an oxygen. This contribution has the ferromagnetic sign and arises as an 
exchange interaction between the hole states (Q) due to their nonlocal nature. The constant 
V is defined in Eq. (|]), h\ is in Appendix. The second term in Eq. ([I]) represents the 
correction to the energy due to the virtual transition of hole from the state Eq. (f|) into 



the two-hole states and back |15||. Here n numerates the two-hole states; matrix elements 



of transitions D n were calculated in Ref. |)15||, AE n are differences in energies between the 



energy of the vacancy and two-hole states at neighbor sites and the ground state energy (see 
Fig. [I]) ; the coefficients x n = 4 for the singlet and x n = —2 for triplet two-hole states. 
The most general expression for the charge-transfer energy is: 

E = E N ~ 1 - E N + E N+1 - E N (8) 

where N refers to the total number of electrons, E^ is the ground state energy, E^^ 1 is the 
minimal energy of a system with one removed and added electron. For our system Eq. (||) 
can be specified as 

E gap = E° gap + AE e + AE h , (9) 

where E® ap is the difference in energies between a singlet and vacancy at local states separated 
by large distance and the ground state (see Fig. ||); AE e and AE h are depths of bands for 
electron and hole (vacancy and singlet). Eg ap can be calculated in the framework of the 
three-band model, whereas for the calculation of AE e and AE^ we will use the t-J model. 
The expression for E® in terms of Eq. (|6|) is very simple: 

E° gap = E c + E v . (10) 

The values of AE e and AEh can be determined from numerous analytical and numerical 



calculations |16| - p9| , [|^,|g,|3)il)il] °f the dispersion relation e(k) for one hole in the t-J 
model on an antiferromagnetic background. There is a general agreement that the hole (or 
vacancy) on the antiferromagnetic background creates a magnetic polaron of a small radius 



j~7| , |2"6|1 , or, in the other words, the carriers are strongly dressed by the spin waves |[23| . The 



influence of antiferromagnetizm and strong correlations are manifested in a special form of 
dispersion relation e(k). For our calculations we use the results from the earlier work by 



Sushkov Ref. [18| where hole wave function and e(k) were obtained variationally: 



e(k) = 1.32J 



1 



- ( AJ - VA2 J 2 + 16t 2 [(l +y)-(x + yhl 



and 

AE = 1.32J + ~ (AJ - ^A 2 J 2 + 16t 2 (l + y)Y (11) 

where for Neel background A = 1.33, a; = 0.56, y = 0.14. Loss of energy due to the broken 
AF bonds (four per carrier) is included. Result for the bottom of the band at t/J = 
2.5 from Eq. ( |Tl~l) AE = —1.2t coincides almost exactly with the recent results of the 



Green function Monte Carlo calculation by Dagotto, Nazarenko, and Boninsegni [23] AE 



— 1.255t. Formula (|TT|) is quite good up to t/J « 5 [18]. 

Let us discuss the characteristic values of all essential parameters which determine the 
observable quantities J and E gap (H), @. In the realistic region of parameters of the three- 
band model one can easily obtain the experimental values J = 0.14 eV and J = 0.17 eV 
for lantan and yttrium systems. The value for E® (|10D was obtained in Ref. [15j and is 



equal to 3.2 eV in the same region of parameters. The characteristic values of AE e and 
AEh fllTD depend on ratios t/J for electrons and holes. These ratios weakly differ and for a 
typical value of t/J = 2.5 we have AE e « AEh ~ 0.42 eV. Thus, the overall gain in energy 
due to magnetic polaron formation is of the order of 1 eV, which is comparable with the 
experimentally observed E gap ~ 2.0 ±0.1 eV. Therefore, the magnetopolaron effect gives 
essential contribution in the value of the charge-transfer gap. 

In the paper p|| it was proposed that the usual phonon polaron effect contributes in the 



observable values of the charge-transfer spectrum. The corresponding gain in energy was 
estimated as 0.5 eV. However, we suppose that due to the Frank-Condon principle the usual 
polaronic effect does not contribute in the optical transition. Magnetopolaron effect has no 
such restriction since it involves electronic degrees of freedom only. We will return to this 
question in Discussions. 

IV. PARAMETERS SENSITIVITY 

Thus, we find superexchange constant J (|7j) and charge-transfer gap E gap ([]) ( [[7|, [2"4"| , pT|| ) 
as functions of the three-band model parameters: 
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" -J \J>pdi "pp., IS, (J d, Up, Vpdji 

Egap = Eg ap {t p di tpp, A, Ud, Up, V p d). (12) 

Both observable quantities are strongly depend on hopping integrals and A = e p — e^, that 
provides the way of fixing these latter by the experimental values of the first. 

As was discussed earlier, abundance of the parameters makes questionable the calculation 
of the effective hopping integral for the t-J model from the three-band model for real CuC>2 
planes. While Coulomb repulsions are known with a fair precision (30% — 50%), the situation 
is complicated due to a very low precision of the direct determination of t pd , t pp and A , 
which mostly affect all effective parameters. Previously, the above parameters have been 
determined from the analysis of spectroscopic data MM . In our recent work we fitted A to 



experimental value of J [15]. 



Now, on the basis of a better understanding of the charge-transfer process and more 
accurate calculations, fixing of the worse known parameters using experimental values of J 
and Eg ap suggests itself. We will show that this procedure keeps effective th inside a narrow 
enough region. 

Firstly, for further discussion we define 0-0 hopping as t pp = jt p d- In order to character- 
ize the above mentioned strong dependence of the parameters ( |12"D one can calculate A(£ p( j) 
at fixed J or at fixed E gap , with other parameters (Coulomb repulsions and 7) as constants 
in the realistic region. We evaluate A vs t p d at J = 140 meV and 170 meV, E gap = 2.0 eV 
and 2.5 eV (see Fig. ^]). Note that the profiles of the curves resemble those at the diagram of 
U(t) for the simple one-band Hubbard model where E gap = U + 2W, W = —at, J = 4t 2 /U, 
and the crossing point uniquely determines U and t. 

To be more specific, we firstly determined A for constant J at an arbitrary t p d and further 
moved up or down along the curve J = const to fix the value of E gap . We used the data 
for La 2 Cu0 4 J = 140 meV J25| and E gap = 2.1 eV (photoconductivity) JEJ. Figures [§§,§ 



show the parameter of our prime interest: the effective integral for hole in the t-J model. 
Parameter 7 = t pp /t p d is 0.5, 0.7, 0.3 for Figs. f|,§,§], respectively. In all figures simple dotted 
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curve corresponds to V P d = U p = 0, dotted ones with crosses V p d — 0, U P = 3 eV, 6 eV, 
dotted ones with triangles V p d = 0.5 eV 1 eV, U p —0, and full curves correspond to including 
both Coulombs V p d = 0.5 eV, U p = 3 eV (upper), Vpd = 1 eV, U p = 6 eV (lower). The 
maximum on the first three curves is due to transition from A > Ud (unrealistic range) to 
A < Ud- All variations of th{Ud, V p d, U p , 7) actually show only weak dependence, and in the 
most preferential region, when all Coulombs are included, th lies between (2.4 -h 2.7) J. We 
believe, that our consideration is quite accurate and well justified. Hence, one can hope that 
the interval for t/J obtained above provides the basis for quantitatively correct calculations 
in the framework of the t-J model. For example, for the recently proposed mechanism of 
superconductivity in the t-J model which provides a very t/J-sensitive (exponentially) gap 
value |2(| . 



Also some other features can be achieved. Figures [7|||§ represent effective hopping \t e \ 
for vacancy, t p d and A respectively, vs other parameters. Here always 7 = 0.5. Strong 
support of our fixing procedure is that selfconsistently determined t p d and A (Figs. |],|9|) he 
in the most appropriate region. From our calculation t p d = 1.2 -f- 1.4 eV, A = 2.5 -j- 4.5 eV 
, that is really close to cluster calculation of Eskes and Sawatzky 0] and to the results of 
other groups p|, ft2[|35 . 



V. EXCITONIC STATE 

The problem of consistently taking into account the Coulomb interaction of the carriers 
in the framework of the t-J - like model remains open. Several recent works are devoted to 



this problem [J38fl . 

As was suggested earlier jj]], the short-range part of this Coulomb interaction may be 
kept by inclusion of the nearest-neighbour Cu-0 repulsion. Since V p d is included in our 
effective model, one can expect that some effects will be caught. 

In the recent work J$J a kink in the optical reflectivity somewhere below the charge- 
transfer peak (at 1.75 eV) was observed. Since it had no associated photoconductivity, it 
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was related to creation of exciton. An essential role of the short-range Coulomb interaction 



was also discussed in Ref. [ZB|]. In our way of reasoning the exciton state, if it exists, is the 
Frenkel exciton, because of short-range nature of interactions. 

Now it is evident that an effective attraction can result from pure magnetopolaron effect 
in the t-J model [|Tj],[J(J. As was shown in works [^,^], "contact" interaction of two holes 
(without charge) dressed by spin fluctuations is attractive for a special symmetry of the 
wave-function. However, the associated energy is very small (< J/3 ~ 0.04 eV). 

It is possible to combine the ideas of the short-range Coulomb and magnetopolaron 
effects. The difference between E® (|10"D when hole and electron are separated and AE n 
(Fig. H) (where n denotes the lowest singlet) when they are close, is the effective Coulomb 
attraction of the "bare" hole and electron (singlet and vacancy). We write an addition to 
Eq. (H) as 

AH C = -V C J2 «3, (13) 

<W>,a/3 

where n e is the electron number operator, n h is the hole number operator. We have found 
that V c is really almost independent from Ud and U p and V c m 0AV p d- "Bare" electron-hole 
attraction itself does not mean Frenkel's exciton effect. One has to show that "dressed" 
electron and hole placed closely possess the lower energy than the mobile ones. In order to 
specify the magnetopolaronic language we reproduce here the wave function of the magnetic 
polaron from the Ref. |18J . For the Ising background it has the more simple form: 



V'fk = —j== Y, rf «T exp(ikr n ), 

yjN/2 n 

dnl = V K\ + l*S+ Y, h n'V ( 14 ) 

n'e<n> 

where n G sub lattice with the spin s = —1/2, h+ a primarily hole operators; at t/J > 1, 
v 2 w 1/2, /i 2 w 1/8. Thus, ansatz consist of mixture of "bare" hole and holes with one 
overturned spin. Contact interaction of these polarons with opposite spins was considered 
in Ref. J5?J (see Fig. |T^). The gain in energy for "attracting polarons" arises from the 
pure t-J model effects |36[ and effective Coulomb attraction (|13"|). The loss in energy is 
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from restriction of the mobility. Competition of these evident effects (without Coulomb 
interaction) provides bound states up to t/J = 2 -h 3 [[!5|,[J1J. Simply acting in a spirit 



of magnetic polaron interaction |36|] we obtain AE exc m 0.35V^. Thus, at V p d = 1 eV, 
AE exc = 0.14 eV. 

More accurate variational construction of the exciton-magnetopolaron wave-function (on 
the Ising background) yields AE exc w 0.5V C (at t/J = 2.5 -S- 3). This wave-function consists 
of the mixture of "bare" hole and electron at neighbor sites and hole and electron with 
overturned spins. It is schematically shown in Fig. |TTla. Thus, at V p d = 1 eV, AE exc ~ 0.2 



eV, that is slightly less than the observed AE e ^P c = 0.25-=-0.35 eV |[28|| . In our calculations the 



interaction between next-nearest-neighbor magnetic polaron was neglected, which produced 



small effect for the pure t-J model [36 but may be essential for the problem with attraction. 



(These configurations are shown in Figs. [TT|b,c). Also, the answer may partly lie in the rest 
of the long-range Coulomb interaction. 

VI. DISCUSSIONS 

The detailed quantitative consideration of some of the effective parameters of the low- 
energy models related to description of the high - T c superconductors presented in this work 
relies heavily on our earlier works. In these works consistent mapping of the three-band 
Hubbard model onto the effective t-J model ( [|l2|,[n||l5[]) has been produced. Taking into 
account all essential interactions enables us to correctly calculate local energies of various set 
of states with different number of particles and matrix elements of interesting transitions. 
Combination of properties of the local "bare" hole and electron (ZR-singlet and vacancy) 
and their magnetic polaron nature as the carriers allows us to approach the calculation of 
some observable quantities adequately. 

We have calculated the superexchange J and charge-transfer gap E gap . Their experimen- 
tal values strongly constrict a possible variation interval for the quantity of great interest: 
the t/J ratio in the t-J model. Selfconsistent calculation of this ratio for a wide range of pa- 
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rameters places it into the region t/J = 2.4-^2.7. Narrowed ranges for the three-band model 
parameters have also been determined: t p d = 1.2 -j- 1.4 eV, A = 2.5 4- 4.5 eV. They coincide 
quite well with earlier cluster calculations, that supports our selfconsistent procedure. An 
excitonic state of the Frenkel type induced by the short-range Coulomb interaction with the 
energy lower than the charge-transfer transition approximately by 0.2 eV is found. 



We have also compared the width of the peak in the 62(0;) at 2.3 eV from Ref. |28] which 
is of the order of 0.5 eV with the total width of the charge-transfer spectrum. This total 
width is equal to the combined width of the vacancy and singlet bands. According to Eq. 
(PI), the width of the hole band is Wh = 2.0 J at th = 2.55 J and width of the electron 
band is W e = 2.2J at t e = 2.15.1 . Resulting total width of the charge-transfer spectrum is 
about 0.6 eV. Thus, the narrowness of the e 2 (uj) spectrum can also be easily reproduced by 
magnetic polaron language. 

One of the essential questions for the Cu02 - planes systems, which we only briefly 
touched, is the phonon polaron effect. If it does not have a projection on optics, or, as 
we believe, the Frank-Condon principle is applicable, our calculated th/ J ratio is the upper 
limit of the real parameter. This is due to mass renormalization for real carriers. There 



is another view on the polaron effect (see Ref. |2l||3"9] ]). It is stated that electron-phonon 
interaction lies in an intermediate range and thus the Frank-Condon principle is not obeyed. 
If such is the situation, we underestimate the depth of the bands (or overestimate E gap ), 
and effective th should be increased. Our estimation shows that this increasing of th is no 
more than 30%. Naturally, this problem requires additional investigations. 
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APPENDIX: 

In this Appendix we present some details of the technical treatment of the problems 
discussed in the paper. According to our previous work [[0J, we use the following transfor- 
mation from the primary oxygen p\ x ,p\ y operators to the operators q\, q~\ of the symmetrical 
and antisymmetrical oxygen states: 

(<7i><Zi) = ^2\pu a cos(k x /2) ±p kb cos(k y /2)] 

k 

x(l + 7 kr 1/2 exp(*kl), (Al) 

where p^ a is the Fourier image of p\ x for q\ and p\ y for qi; p k6 is the Fourier image of p\ y for 
q\ and p\ x for q\\ 7 k =(cos(k c a)+cos(k, / a))/2. The summation in Eq. ( |A1| ) is produced over 
the Brillouin zone, and the lattice constant a= 1. 

Since the groundstates of both undoped and doped systems do not consist of antisymmet- 
rical oxygen state [P,[F|,|i"5] , the reformulated Hamiltonian (0,0) where only essential degrees 
of freedom are kept, is conveniently expressed through the local and hopping parts [I5|. The 
local part is: 

Hioc = e d Y, nf a + (e p - fi t pp ) £ nf a 

l, a l,a 

+ U d Y^ « + V pd f £ n d la n% 

l l,a/3 

+ U p h J2 n h n li + 2t P dA ^2{dt a qia + H.c.) , 

l l,a 

AH int = Vpdfi J2 n t n h 

<W>,af3 

- 2U p h Y, (SfSf, - -nlnl) , (A2) 

<w> 4 



1 
2* 



with Sf = ^qt a cr a pqip, n q = n q n + n\ v (A3) 
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The hopping part is: 



H hop = 2t pd Xi J2 ( d ia<li'a + H - c -) 

<W>,a 

- Ztppflx 22 Qla^l'a, 
<//'>, a 

AH hop = V pd f J2 nf a lq+ qi , p + H.c.} 

<ll'>,a/3 

+ U p h' £ <[ g +^ + F.c.]. (A4) 

<W>,a 



All constants X,fi,f,h in Eqs. ( [A2| , [A4] ) are of Wannier nature. Their Fourier images and 



magnitudes are given in Ref. fL5| |. In order to group them together we reproduce 



A =0.9581 Ai=0.1401 /i =1.4567 /ii=0.2678 
/o=0.9180 /i=0.2430 /i =0.211 /ii=0.059 
/' =0.1342 /i'=0.030 . 

We have treated the Hamiltonian H\ oc + AHi nt ( |A2| ) in the selfconsistent mean- field approx- 
imation (JT5J that enables us to solve the problem of the local states at site with different 
number of holes. The matrix elements of the Hamiltonian ( |A4l ) between states with a singlet 
or vacancy at different sites in initial and final states lead to the following expression for the 
hopping constants: 

t h = 2t\ 1 (W 1 V - \f2U x U')(yi x XJ' - VWxV) 
+ t p fi 1 (W 1 U' -V2V 1 V') 2 /2 

- (ynftwji' - VzUphtVtV'KWxU' - VwiV), 

t e = - UX^'V" - tp^iV") 2 , (A5) 

where U', V are the coefficients of |/ > - state nearest to singlet, U", V" - those nearest to 
vacancy. The coefficients U', V and U", V" are slightly different from the ones in Eq. (|J) 
due to a distortion of \f > - states by the nearest vacancy or singlet. This distortion has its 
origin in the short-range Coulomb repulsions due to U p and V p d terms in the Hamiltonian 
(@). Note that it changes the energy of the \f > - states nearest to a vacancy or singlet. We 
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take this effect into consideration when we calculate the quantity E® . Thus, we have for 



the total energies E c and E v for a singlet and a vacancy Eq. ( [10] ) 



E c = E s + AE fs - 5E f , 

E V = E + AE f0 - 5E f , (A6) 

where E s and Eq are the local energies of singlet and vacancy; Ef 8 and Efo are the energies 
of the |/ > - states nearest to a singlet or vacancy; Ef is the energy of |/ > - state in an 
undoped sample. These states are schematically shown in Fig. ||. 

When the vacancy and the two-hole state are created at the neighbour sites (see Fig. 
[I]) the difference in energy between this state and the ground state is determined by the 
relation 

AE n = E n + E Q + 3E fs + 3E f0 - 8E f . (A7) 

These energies AE n are involved in calculation of the superexchange constant J (Eq. (0)) 
and in the energy of the Coulomb attraction of a singlet and vacancy in Sec. V. 
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FIGURES 
FIG. 1. The nearest-neighbor two-hole state and vacancy. Black circle denotes two-hole state, 

empty circle denotes vacancy. Crosses are one hole states (spins). 

FIG. 2. The separated ZR singlet (hole) and vacancy (electron). Black circle denotes singlet, 
empty circle denotes vacancy. Crosses are one- hole states (spins). 

FIG. 3. A vs t p d at constant J or E gap . Ud = 7 eV, U p = 3 eV, V p d = 1 eV, t pp /t p( i = 0.5. 
Full curve - J = 140 meV, dashed curve - J = 170 meV, full curve with markers - E gap = 2.0 eV, 
dotted curve - E gap = 2.05 eV. 

FIG. 4. Effective hopping integral for t-J model hole vs Ud- Dotted line - V p d = U p = 0, dotted 
with crosses - V p d = 0, U p = 3, 6 eV, dotted with triangles - U p = 0, V p d = 0.5, 1 eV, full curves - 
V p d = 0.5 eV, Up = 3 eV (upper), V p d = 1 eV, U p = 6 eV (lower), 7 = 0.5. 

FIG. 5. All notations as for Fig.4; 7 = 0.7. 

FIG. 6. All notations as for Fig.4; 7 = 0.3. 

FIG. 7. Effective hopping integral for electron in the t-J model vs Ud, curves markers as for 
Fig.4; 7 = 0.5. 

FIG. 8. Cu-0 hopping integral vs Ud, curves markers as for Fig.4; 7 = 0.5. 

FIG. 9. A vs U d ; 7 = 0.5. 

FIG. 10. Configuration of interacting magnetic polarons from Ref. [36]. 

FIG. 11. (a) Exciton-magnetopolaron; (b),(c) next-nearest-neighbor magnetic polarons. 
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